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We have measured the process B ± -> (K* ± -> ^tt )^ with 232 million T(4S) -> BB decays 
collected with the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. From a 
signal yield of 89 ± 26 events we obtain the branching fraction B(B ± — > K* ± tt°) = [6.9 ± 2.0(stat) ± 



4 



1.3(syst)] x 10 6 with a statistical significance of 3.6 standard deviations including systematic 
uncertainties, and a charge asymmetry of 0.04 ± 0.29(stat) ± 0.05(syst). 

PACS numbers: 11.30.Er, 13.25.Hw 



Branching fraction and CP- asymmetry measurements 
of charmless £?-meson decays provide valuable con- 
straints for the determination of the unitarity triangle 
constructed from elements of the Cabibbo-Kobayashi- 
Maskawa quark-mixing matrix [HQ. They test the accu- 
racy of theoretical models such as those based on QCD 
factorization Q or SU(3) flavor symmetry £|. It has 
been argued that the influence of final-state interactions 
like charming penguins |5j,|6j,|jJ and similar long-distance 
rescattering effects Q on both the branching fraction and 
CP asymmetry of B — Kn decays may be significant. 
In this context, the decay B^ — > K* ± ir a is particularly 
interesting in the light of recent measurements of direct 
CP-violation in the 5° ->■ and B ± -> K ± tt° chan- 

nels [1MII3- Comparison to the B ± -> if*°7r ± decay 
mode |l2j can provide information about the dominance 
of penguin diagrams. Here we present a measurement 
of the branching fraction B(B ± -> if *(892) ± tt°) and its 
charge asymmetry 

_ N(B- -> K*-ir°) - N(B+ -> K*+ir°) 
CP ~ N{B- ^ K*-tt°) + N{B+ K*+tt°) 

based exclusively on B ± decays to the K^tt^-k final 
state. The data used in this analysis were collected 
with the BABAR detector [l3| at the PEP-II asymmetric- 
energy e + e~ storage ring at SLAC. Charged-particle 
trajectories are measured by a five-layer double-sided 
silicon vertex tracker and a 40-layer drift chamber lo- 
cated within a 1.5-T solenoidal magnetic field. Charged 
hadrons are identified by combining energy-loss informa- 
tion from tracking with the measurements from a ring- 
imaging Cherenkov detector. Photons are detected by a 
CsI(Tl) crystal electromagnetic calorimeter with an en- 
ergy resolution of a E /E = 0.023(^/GeV)^ 1 / 4 © 0.014. 
The magnet's flux return is instrumented for muon and 
K® identification. 

The data sample includes 232 ± 3 million BB pairs 
collected at the T(4S) resonance, corresponding to an 
integrated luminosity of 211 fb _1 . It is assumed that 
neutral and charged B meson pairs are produced in equal 
numbers In addition, 22 fb _1 of data collected at 

40 MeV below the T(AS) resonance mass were used for 
background studies. We performed full detector Monte 
Carlo (MC) simulations equivalent to 460 fb _1 of generic 
BB decays and 140 fb _1 of continuum quark-antiquark 
production events. In addition, we simulated over 30 
exclusive charmless B decay modes, including 1.2 million 
signal B ± — > K* ± tt° decays. 
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B meson candidates are reconstructed from one 
charged track and two neutral pions. The charged 
track used to form the _B ± — * K*^^ candidate is re- 
quired to have at least 12 hits in the drift chamber, to 
have a transverse momentum greater than 0.1 GeV/c, 
and to be consistent with originating from a £>-meson 
decay. Its signal in the tracking and Cherenkov detec- 
tors is required to be consistent with that of a kaon. The 
kaon selection algorithm is 70-92% efficient within the 
relevant momentum range, with a misidcntification rate 
of less than 7%. We remove tracks that pass electron 
selection criteria based on dE/dx and calorimeter infor- 
mation. Neutral pion candidates are formed from two 
photons, each with a minimum energy of 0.03 GeV and 
a lateral moment [l5| of their shower energy deposition 
greater than zero and less than 0.6. The angular ac- 
ceptance of photons is restricted to exclude parts of the 
calorimeter where showers are not fully contained. We 
require the photon clusters forming the tt° to be sepa- 
rated in space, with a ir° energy of at least 0.2 GeV and 
an invariant mass between 0.10 and 0.16 GeV/c 2 . 

Two kinematic variables, AE = E* B — y/s/2 
and the beam energy substituted mass tties = 
\J (s/2 + po • Pb) 2 /E'q — Pg, ar e used for the final se- 
lection of events. Here Eg is the i?-meson-candidatc 
energy in the center-of-mass frame, Eq and yfs are the 
total energies of the e + e~ system in the laboratory and 
center-of-mass frames, respectively, and po and ps are 
the three-momenta of the e + e~ system and the B candi- 
date in the laboratory frame. For correctly reconstructed 
K* ±r K® candidates AE peaks at zero, while final states 
with a charged pion, such as J5 ± — > p ± 7r°, shift AE by 
approximately 80 MeV on average. Events are selected 
with 5.20 < m ES < 5.29 GeV/c 2 and \AE\ < 0.20 GeV. 
The AE limits help remove background from two- and 
four-body B decays at a small cost to signal efficiency. 

Continuum quark-antiquark production is the domi- 
nant background. To suppress it, we select only those 
events where the angle #gph m ^ ne center-of-mass frame 
between the direction of the i?-mcson-candidatc and 
the sphericity axis of the rest of the event satisfies 
| cos 6*g ph | < 0.9. In addition, we construct a non-linear 
discriminant, implemented as an artificial neural network 
(ANN) that uses three input parameters: the zeroth- and 
second-order Legendre event shape polynomials Lq, Li of 
the momenta and polar angles of all candidates in the 
rest of the event, and the output of a multivariate, non- 
linear B-meson-candidate tagging algorithm [l6j. ANN 
is peaked at 0.5 for continuum- like events and at 1.0 for B 
decays. We require ANN > 0.58 for our event selection. 
To further improve the signal-to-background ratio we re- 
strict the effective invariant mass of the K* candidate to 
0.8 < m K7r < 1.0 GeV/c 2 . Neutral -pion combinatorics 
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lead to 30% of our signal events having more than one 
candidate per event. We choose the best candidate based 
on a x 2 formed from the measured masses of the two 7r° 
candidates within the event compared to the known ir° 
mass |l7|. 

After the selection described above, the B ± — > K*^ir° 
selection efficiency is 16.5%. In MC studies, the signal 
candidate is correctly reconstructed (64.5 ± 6.5)% of the 
time. The remaining candidates come from self-cross- 
feed (SCF) events, which stem primarily from swapping 
the low energy ir° from the resonance with another from 
the rest of the event. The fraction of SCF events in which 
the track was swapped with an oppositely charged track 
was found to be negligible. 

MC events are used to study backgrounds from other 
B-meson decays. The dominant contribution comes from 
b — > c transitions; the next most important is from 
charmless _B-meson decays. The latter tend to be more 
problematic as the branching fractions are often poorly 
known, and because they may peak at the same invariant 
mass as the signal B ± — ► K*^7r° events. Thirteen indi- 
vidual charmless modes show a significant contribution 
once the event selection has been applied. These modes 
are added into the fit fixed at the yield and asymmetry 
determined by the simulation. Wherever branching frac- 
tions are not available, we use half the upper limit. If no 
charge asymmetry measurement is available, we assume 
zero asymmetry. 

Although all other known resonant K* states - sub- 
sequently referred to as K** - lie outside our if* (892) 
mass window, some may still contribute due to their large 
width. To estimate the contribution to the signal we se- 
lect a region in the K ± n° invariant mass between 1.2 
and 1.6 GeV/c 2 , motivated by the presence of the broad 
Kq(143Q) resonance which decays predominantly to Kir. 
In this region we make a full maximum likelihood fit to 
the three variables AE, mEs and ANN in an analogous 
way to how we fit our signal (see below), and extrapolate 
the result of this fit to the if* (892) signal region using a 
B ± -> #£(1430)%° MC. The fit to the 2^(1430) region 
yields 263 ± 34 events, which translates to 34 K** events 
contributing to the background in our signal region. We 
assign a 100% systematic uncertainty to this number to 
cover possible interference effects as well as uncertainty 
arising from the lineshapes of K** resonances, which are 
not well established. 

The non-resonant B ± — > K^tt ^ branching fraction 
has, to date, not been measured. To estimate the signif- 
icance of its contribution we select a region of the Dalitz 
plot of rn? K7T - for the ir° from primary and secondary de- 
cay - that is far from the signal as well as if* (1430) and 
higher K** resonances and which has low levels of con- 
tinuum background. A likelihood fit in this region yields 
6±8 events, which translates into less than three events in 
our K* signal region, assuming the non-resonant events 
are distributed evenly across the Dalitz plot. We conse- 
quently deem the non-resonant contribution negligible. 

An unbinned maximum likelihood fit to the variables 



tties, AE, rriK-K, and ANN is used to extract the total 
number of signal B^ — > K*^7T and continuum back- 
ground events and their respective charge asymmetries. 
The likelihood for the selected sample is given by the 
product of the probability density functions (PDF) for 
each individual candidate, multiplied by the Poisson fac- 
tor: 

1 N 
C=^e~ N ' [NTH*, 

i=l 

where N and AT' are the number of observed and expected 
events, respectively. The PDF Vi for a given event i is 
the sum of the signal and background terms: 

+ (l-Q t A Si z)(l-f)V? is ] 

+E< kSx ^ 1 -^ Bks )^ ks ' 

3 

where Qi is the charge of the kaon in the event, 
jySig^Bkg-j and ^Sig^Bkg^ are the yield and asymmc _ 

try for signal and background component j, respectively, 
and / = 35.5% is the fraction of SCF signal events. 
The j individual background terms comprise continuum, 
b — > c decays, K**, and 13 exclusive charmless B de- 
cay modes. The PDF for each component, in turn, is 
the product of the PDFs for each of the fit input vari- 
ables, V = 'P mESt AE'PANN'P mKn . Due to correlations be- 
tween AE and tties, the V mESi AE for signal and all back- 
ground from B decays are described by two-dimensional 
non-parametric PDFs obtained from MC. For con- 
tinuum background, the correlations in AE and tties are 
~ 1%, hence a separate PDF is used for each of them; 
?ties is well described by an empirical phase-space thresh- 
old function [l9| and AE is parameterized with a sec- 
ond degree polynomial. The parameters of the contin- 
uum PDFs are floated in the fit. ANN is described by 
a non-parametric PDF for continuum background and 
by a Crystal Ball function |20| for all other modes. For 
y mK T, i one-dimensional non-parametric PDFs obtained 
from MC are used to describe all modes except the sig- 
nal mode itself, which is described by a Breit-Wigner 
line-shape combined with a first degree polynomial. The 
parameters for this PDF are held fixed to the MC values 
and varied within errors to estimate systematic uncer- 
tainties. 

A number of cross checks confirm that the fit is un- 
biased. In 1000 separate MC experiments we generate 
the expected number of events for the various fit compo- 
nents before using the maximum likelihood fit to extract 
the yields and asymmetries. The distributions for each 
component are generated from the component's PDF, 
giving values for the fit variables mss, AE, ANN, and 
rriK-K- The expected number of events is calculated from 
the branching fraction and efficiency for each individual 
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TABLE I: Breakdown of systematic uncertainties. 
Absolute systematic uncertainties on yields 



Source 


a™ d (Events) 


Background normalization 


±14 


PDF shapes 


+ 2.1 
- 4.0 


SCF fraction 


± 1.8 


AE shift 


± 2.2 


Total 


±14 



Relative systematic uncertainties on B(B ± 
Source o"s ys t.(%) 



K TT 



Efficiency estimation 
B counting 



±7.3 
±1.1 



Total 



±7.4 



Systematic uncertainties on Acp 



Source 



A CP 

J Syst. 



Background normalization 
Detector asymmetry 
Background asymmetry 



+0.018 
-0.010 

±0.003 

+0.049 
-0.041 



Total 



+0.054 
-0.043 



mode. The generated number of events for each fit com- 
ponent is determined by fluctuating the expected num- 
ber according to a Poisson distribution. The test is re- 
peated using samples with differing asymmetry values. 
We repeat these MC studies using fully simulated sig- 
nal B^ — ► K*^^ events instead of generating the sig- 
nal component from our PDFs. This verifies that the 
signal component is correctly modeled including correla- 
tions between the fit variables. Finally, omitting itlktv 
as a fit variable has no significant influence on the signal 
yield, indicating that our treatment of K** background 
is indeed effective. 

Individual contributions to the systematic uncertainty 
are summarized in Tabled We calculate the uncertainty 
of the continuum background estimation directly from 
the fit to data. The backgrounds from B decays are de- 
termined from simulation and fixed according to their 
efficiencies and branching fractions. For those individ- 
ual decay modes which have been measured we vary the 
number of events in the fit by their measured uncertainty. 
For all others we vary the amount included in the fit by 
±100%. For the b — > c component we fix the rate based 
on the number calculated from MC samples and vary the 
amount based on the statistical uncertainty of this num- 
ber (6%). The shifts in the fitted yields are calculated 
for each mode in turn and then added in quadrature to 
find the total systematic effect. The largest individual 
contribution comes from the K** estimation. 

To take into account the variation of the two- 
dimensional non-parametric PDFs used for AE and ttieSj 
we smoothen the MC-generated distributions from which 
the PDFs are derived. For tuk-k and ANN , the parame- 
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FIG. 1: Likelihood projection plots for the four fit variables, 
(a) mES, (b) AE, (c) ANN, and (d) mien- In each plot the 
solid line represents the total PDF, the dotted line represents 
the total background, the dotted-dashed line represents the 
continuum contribution, and the dashed line represents the 
signal component. The plots contain a subset of the events 
defined by a likelihood ratio of at least 0.1 (see text). 



terizations determined from fits to MC events are varied 
by one standard deviation. The systematic uncertainties 
are determined using the altered PDFs and fitting to the 
final data sample. The overall shifts in the central value 
are taken as the size of the systematic uncertainty. 

We vary the SCF fraction by a conservative estimate 
of its relative uncertainty (±10%) and assign the shift 
in the fitted number of signal events as the systematic 
uncertainty of the SCF fraction. 

To account for differences in the neutral-particle re- 
construction between data and MC simulation the signal 
PDF distribution in AE is offset by ±5 MeV and the data 
refitted. The larger of the two shifts in the central value 
of the yield is 2.2 events, which is taken as the systematic 
uncertainty for this effect. 

Corrections to the 7r° energy distribution, determined 
using various control samples, add a systematic uncer- 
tainty of 7.2%. A relative systematic uncertainty of 1% is 
assumed for the kaon identification. A relative systematic 
uncertainty of 0.8% on the efficiency for a single charged 
track is applied. Adding all the above contributions in 
quadrature gives a relative systematic uncertainty on B 
of 7.3%. Another contribution of 1.1% comes from the 
uncertainty on the total number of B events. 

The cross section for the interaction of kaons with pro- 
tons and neutrons differs with charge. At low momenta 
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this can introduce a bias to the observed charge asymme- 
try. We estimate this bias by modelling the average loss 
of kaons from a sample based on the K* ± n° signal MC 
using the known detector material constants, and find 
Akp = —0.0031 ± 0.0006, which is negligible compared 
to the precision at which we measure Acp- 

To calculate the effects of systematic shifts in the 
charge asymmetries of background modes, each mode is 
varied by its measured uncertainty. For contributions 
with no measurement, we assume zero asymmetry and 
assign an uncertainty of 20%, motivated by the largest 
charge asymmetry measured in any mode so far [t| . The 
individual shifts are then added in quadrature to find 
the total systematic uncertainty. The greatest individual 
contribution comes from the K** estimate. In addition, 
the effect of altering the normalizations of the B back- 
grounds affects the fitted asymmetry. The size of the shift 
on the fitted Acp is taken as the size of the systematic 
uncertainty. 

A total of 23,465 events were fitted, of which 11,960 
had positively charged candidates. The central value of 
the signal yield from the maximum likelihood fit is 89 ±26 
events, over an expected background of 634 ± 40 events 
from other B decays. We obtain a branching fraction of 
B(B ± -> (K^ -> i^Tr )?!- ) = [2.31 ± 0.67 ± 0.42] x 
10 -6 and charge asymmetry of AcpiB^ — > (K*^ — > 
i4r ± 7r°)7r°) = 0.04 ± 0.29 ± 0.05, where the first error 
is statistical and the second one systematic. Alterna- 
tively, we calculate the 90% confidence upper limit on 
the B ± -> (iT* ± -> i^Tr )?!- branching fraction to be 
3.9 x 10 -6 . Compared against the null hypothesis, the 
statistical significance y— 2 ln(£ Null / Cm ax) of the yield 
amounts to 4.1 standard deviations. The fit was redone 
fixing the signal yield to the lowest yield allowed account- 
ing for all possible combinations of systematic uncertain- 
ties. The significance of this result corresponds to 3.6 
standard deviations. 



The results of the fit are illustrated in Fig. ^ The 
plots are enhanced in signal by selecting only those events 
which exceed a threshold of 0.1 for the likelihood ra- 
tio R = (jVSigpSigy^SigpSig +£ < JVf k8 Pf ks ), where 
N are the central values of the yields from the fit and 
V are the PDFs with the projected variable integrated 
out. This threshold is optimized by maximizing the ra- 
tio S = (N Sis e sig )/(^V Sig e Sig + N? kg ef kg ) where 
e are the efficiencies after the threshold is applied. The 
PDF components are then scaled by the appropriate e. 

In conclusion, we have measured the charge asymme- 
try and branching fraction for the decay B^ — > (K*^ — > 
K^tt )^ using a maximum likelihood fit. Assuming 
a secondary branching fraction of 1/3 for the K*^ — > 
K ± tt° final state our result implies B(B ± — ► K* ± ir°) — 
[6.9 ± 2.0 ± 1.3] x 10~ 6 , and a charge asymmetry of 
Acp = 0.04 ± 0.29 ± 0.05 where the first error is sta- 
tistical and the second error systematic. The statisti- 
cal significance of the branching fraction result including 
systematic uncertainties is calculated to be 3.6 standard 
deviations, showing evidence for this decay. The system- 
atic error of the branching fraction and asymmetry is 
dominated by the contribution of K** resonances. 
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